Shock waves in one-dimensional Heisenberg ferromagnets. 
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We use the SU(2) coherent state path integral formulation with the stationary 



[ — ■ phase approximation to investigate, both analytically and numerically, the existence 

of shock waves in the one-dimensional Heisenberg ferromagnets with anisotropic ex- 
change interaction. As a result we show the existence of shock waves of two types, 
O ■ "bright" and "dark", which can be interpreted as moving magnetic domains. 
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The Heisenberg model is certainly one of the most important models of condensed 
matter physics on which a large amount of work has been done HJ. In spite of this, 
the study of its properties, both at the quantum and at the classical level, is still a non 
exhausted subject of ever-continuing interest. In the one dimensional case the model 
is exactly solvable and it is known that the excitation spectrum consists of quantum 
solitons which in the semiclassical limit can be seen as bound states of many magnons. 
On the other hand, it has recently been shown that nonlinear lattices, besides solitons, 
may support other kinds of excitation which behave at the initial stages of their evolution 
like shock waves in liquids or gases. Such waves have been reported both in integrable 
0-^] and nonintegrable |5|-|| lattices. 

The aim of the present paper is to show the existence of shock waves in the one- 
dimensional Heisenberg ferromagnet. To this end we use the SU(2) coherent state |J 
path integral formulation with the stationary phase approximation JlIJ , to derive classical 
equation of motion out of the original quantum spin model. Within the framework of this 
approximation the system is shown to be described by a discrete nonlinear Schrodinger 
(DNLS) like equation which has a hamiltonian structure and a non standard Poisson 
bracket. Moreover, the classical system preserves the conservation of the analogue of the 
quantum total spin (for isotropic case) operator as well as of its z-projection. We find that 
for suitable condition the excitations of this system naturally display shock waves with 
sharp rectangular profiles moving on uniform backgrounds. Such waves can exist both 
above background (bright shock) and below background (dark shock) and on the contrary 
of other excitations, which may decay into soliton trains or background radiation, they 
are very stable. We remark that similar solutions were found also in a deformable DNLS 
system P,LT| and in a chain of two level atoms describing the propagation of Frenkel 
excitons |T2[ . We give a numerical and an analytical description of these phenomena both 
in terms of dispersion relations and in terms of a small amplitude multiscale expansion. 

The quantum Heisenberg Hamiltonian is written as 
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H = -J2 J(n,m)[(S:S x m + S^J + \S z n S z m ] (1) 

(m,n) 

where S n = S%, S z ) are spin operators of spin magnitude S, J(n,m) is the exchange 
interaction constant and A is the anisotropy of the exchange XF-like (A < 1) and Ising-like 
(A > 1) interactions. In what follows we consider the case of nearest neighbor interaction 
J(n, m) = J ■ (5 n>m+ i + 5 n ,m-i) with J > and we denote, as usual, with = ± iS% 
the raising and lowering operators. To derive classical equation of motion it is suitable to 
use SU(2) coherent states 

in terms of which we write the state of the system ([I]) in the form 

iA)=niM. (3) 

n 

Here complex variable and \ i) n denotes the ground state of a single site. The 

time evolution operator of the Heisenberg hamiltonian between the initial state |A)j at 
time £j and a final state |A) f at time tf can be written in terms of path integral as 

J (A|exp(-^)|A) / = / 3(A)«p(i£) (4) 

where S = /// Ldt with L given by 



Using the stationary phase approximation 5S = in @ one readily obtains the following 
equation of motion 



dt \ l + |/i n+1 | 2 1 + 
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Mn(l — l^n+ll ) , Mn(l — l/^n-ll ) 

-1 



+aj , ^"7' y + , |A ""7' > . (6) 



It is of interest to note that and its complex conjugated follow from Hamilton's 
equations 



with the noncanonical Poisson bracket 

df dg df dg 



{f,9} = -^ + \»n\ 2 T 



_dp n dp n dp n dp r 
and with the classical hamiltonian H c given by (A\H\A) i.e. 

tj 7\~^ I P"n+ll J >n + A (1 — |/^ n | 2 ) (1 — 1/ ' n + J ' \ , n , 



(1- 
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One readily checks that the conservation of the z-component of the quantum total spin is 
reflected in the classical system (|6]) in the conservation of the quantity 

g ^(AlE^|A) = ^E [l:Kv (io) 

In the isotropic case (A = 1) besides ( 110 ) there is another conserved quantity which is 
just the analogue of the total spin of the quantum system 

s2s<a|s2|a) ^5(TTO™ < u > 

[in an infinite chain with |/i n | — > p as n — > ±00, expressions ( |10D and (|Il]) should be 
properly normalized to be finite quantities]. The classical system keeps, therefore, 
some important symmetry of the original quantum model [note also that the dynamical 
variables (p n ,P>n) can be related by an inverse stereographic projection to the motion of 
vectors on a sphere (classical spins)]. In order to investigate shock solutions of (^) it 
is suitable to consider solution propagating against a nonzero background of the form 
pffl = pexp(—iujt + ikn) where p < 1 and 



1 -P 2 

u = J(A-2cos(fc)) hr. (12) 

1 + p 2 



(here and in the following we fix h — 1). The stability of the background can be studied 
with help of the substitution p n = (1 + ip n )pexp(—iut + ikn), where \ip n \ -C \p n \, in 
(^). By linearization we obtain the dispersion relation Q(K) associated with the linear 
equation for ip n [ex exp(— iflt + iKn)] 
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{cos 2 (A;)(l + p 2 ) 2 — cos 2 (fc) cos(if)(l - p 2 ) 2 - 4Ap 2 cos(A;) cos(K)} . (13) 

Thus the background is stable (i.e. f2 is real at all K), if cos k > X and > cos k > — A. 
Naturally, in what follows the analysis will be restricted to this region of the parameters. 
In order to get the equation governing the initial stages of the evolution of a shock wave 
we use the small amplitude expansion \x n = (p + a n ) exp[z(— ut + kn — <p n )], where the two 
real quantities a n and <p n are considered depending on slow variables X = 771, T = jt, 
and t = 7 3 t, (with 7 C 1) and are represented in the form a n = 7 2 a^°- ) + 7 4 o^ + 
4>n — 70l + T 3 0n + ••• • Collecting all the terms of the same order in 7 we arrive at 



a series of equations. In the zero order we recover the dispersion relation (jT2|) . In the 
second and third orders we get the following equations 



8pJ (cosA;- A)a (0) -2Jsin(fc)- — ~^r> (14) 



dT (1 + p 2 ) 2 " ' " J l + p 2 dX 

dai0) T fjr^^ 9 • <1,M 1 -P 2da{0) 

-gjr = PJ™m-^ ~ 2^)7——. (15) 

It is suitable to introduce new variables (£±, T) instead of (X, T), where £± = X — c±T 
and the velocities c± are given by 

2J 



1 — p ) sin k ± pJ 2 cos A;(cos k — A) 



(16) 



(1 + P 2 ) 

Comparing this result with ( JTBj ) one sees that c± = dfl±/dK at = 0, i.e. c± are group 
velocities of two branches of the spectrum in the center of the BZ. Then it follows from 
QJJ), (|TJ) that solutions a® = a (0) (£±) = a± and = (o) (£±) are related by 



VcosA; c^ - 1 



a± = T(l + P 2 ) , Z c ■ (17) 

2J2(cosfc- A) ^± 



For brevity we drop out the explicit form of the equations appearing in the forth and fifth 
orders of 7 and present simply the condition of their compatibility. This condition has 
the form of a Korteweg - de Vries (KdV) equation 



da± da± , M <9 3 a± 

— — + a(k)a±— - + I3{k)^r- = 
or ot,± 



da 



with 



a{k) 



4J 



2\2 



[— 5p sin k + 2A tan k 



[1 + p 



±(1 - 3p 2 W2cosfc(cosfc - A) ± 



AV2 cos 3 /c , 9n 

- 1 " ' xi - L 9 ; , X (3- P 2 ) 

2V cos /e — A 



(19) 



and 



J 



4p(l + p 2 )VcosT 
"2 



r 8 

— — |-psin k\J cos — A(l — p 2 ) 



±V2 cos k 



p (cos — A) — 4p A — cos k(l — p 



2\2 



(20) 



From (gD,(gg) i1] follows that if 



Ay/2 P (1 - p 2 ) sin k\ C ° Sk , A = ±3[(1 - p 2 ) 2 cos + 4Ap 2 ] T 2p 2 (cos jfe - A) (21) 

V COSK 

is satisfied, the coefficient (3(k) becomes zero and the KdV equation reduces to the well- 
known equation 



da± da± 
+ a(k)a±- 







(22) 



which support shock solutions. This implies that for parameter values satisfying fl2"T| ) 
shock wave should develop in the spin chain. 



There are two facts concerning (|22D to be mentioned here. Note that from expression 
(f22|) it follows that at k = there exists only one background < p < 1 at which /3(k) 
equals to zero while this is not true for k ^ 0. Moreover, equation (|22"D is not satisfied for 
all A values but there exists a maximal value \ m ax = 1/7 above which (^) does not have 
physical meaning. 

To check the above predictions, we have numerically integrated on a long chain 
(we neglect boundary conditions), taking as initial condition a bell shaped bright or dark 
pulse of the type 
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(note that with this initial condition rectangular shock profiles should develop). In Fig. 1 
we have reported the profile which develop out from an initial bright pulse of amplitude 
\A\ = 3.6, after an evolution time of T = 420. The background is moving in-phase (k = 0) 
with p = 0.8 and for parameter values given by J = 0.8, and A derived from (|2"T|). From 
this figure we see the appearance of a leading rectangular shock profile followed by solitons 
and background radiation. The shock wave connects the uniform background field with a 
local plateau with two sharp transitions at the edges. If we define the local magnetization 
as M n = (1 — \p n \ 2 ) I (1 + \p n \ 2 ) we have that the local magnetization in the rectangular 
shock waves of Fig.l does not change in time and is different from the surroundings. This 
suggest the interpretation of such solutions as propagating magnetic domains. 

A similar result can be obtained starting from an initial dark profile as shown in 
Fig. 2. In this case the shock plateau develops below the background and therefore it can 
be referred to as a dark shock. Notice that in the above context dark and bright pulses 
correspond respectively to domains with higher and lower magnetization compared with 
the magnetization of the background. 

Following the time evolution of the shock profiles in Fig.s 1,2 we find that the rectan- 
gular waves separate from the other components (solitons and radiation) and stay stable 
over long time. We have numerically checked that (^l|) is a necessary condition for creation 
of shock waves. In Fig. 3 we have reported the evolution profile for the same parameter 
values of Fig.l except for A = 0.185 not satisfying relation (^l|). We see that the rectan- 
gular shock is destroyed and oscillations develop on the wavefront. The same is observed 
for other choices of parameters for which Eq. ( f21~|) is not satisfied. We also checked that 
shocks remain stable upon collision with other excitations (the same was found in ref. 
H) this suggesting in them the presence of a strong soliton component. This leads to the 
interpretation of shocks as bound states of many solitons (in the quantum system they 
should correspond to long "string" excitations i.e. to bound states with a large number of 



quasiparticles). To check this interpretation, however, further investigations are required. 

In conclusion, we have shown both analytically and numerically the existence of shock 
waves in ID Heisenberg ferromagnets. Physically, this is a novel example of manifestation 
of classical fluid dynamics in quantum magnetic systems. Previously, classical behaviors 
in magnetic systems were reported on multi-magnon instabilities and chaos in pure and 
doped YIG and some antiferromagnets [13|. In order to observe such shock waves as 
those obtained here in magnetic systems, generation of macroscopic number of magnons 
is generally required. One of promising candidates to realize this may be highly pummped 
YIG. 
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FIGURES 

FIG. 1. Evolution of a bright shock against nonzero background with k = p = 0.8, and 
for parameter values J = 0.8 and A determined from (|2T|). 

FIG. 2. Same as in Fig.l but for a dark initial condition and for parameter values p = 0.8 
J = 1 and A determined from (pip. 

FIG. 3. Evolution profile for the same parameter values as in Fig.l but with A = 0.185 not 
satisfying Eq. (pl|) . 
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